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This study investigates the spin susceptibihty in the superconducting state of a multilayer system with a layer depen- 
dent Rashba spin-orbit coupling, representing a locally non-centrosymmetric superconductor It is shown that the charge 
imbalance between the different layers yields a strong layer dependence in the susceptibility, most significantly for weak 
spin-orbit coupling, which is in contrast to the situation with equally distributed charge. These results can be relevant 
for experimental test in multilayer high-r,; cuprates as well as in superconducting CeCoIns/YbCoIns superlattices with 
more than two layers. 
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Since the discovery of superconductivity in the heavy- 
fermion compound without inversion symmetry, CePt3Si,'^ 
the non-centrosymmetric superconductivity has been 
studied extensively.^' Subsequently, several new non- 
centrosymmetric superconductors with unique properties 
have identified among heavy-fermion materials'*^' and oth- 
ers. ^"'^' Many exotic properties induced by the antisymmetric 
spin-orbit coupling have been studied in various contexts, 
such as the exotic superconductivity and spintronics.^* 

On the other hand, many other superconductors lack local 
inversion symmetry even while having the global inversion 
symmetry in the crystal structure. In contrast to the uniform 
antisymmetric spin-orbit coupling in non-centrosymmetric 
systems, the antisymmetric spin-orbit coupling appears in 
such systems with special unit cell structures locally loosing 
inversion symmetry. An interesting class of such "locally non- 
centrosymmetric system" is the multilayer systems in which 
the Rashba-type spin-orbit coupling is induced by the local 
violation of mirror symmetry.'^' For instance, the recent ex- 
periment discovered superconductivity in the artificial super- 
lattices of heavy-fermion superconductor CeCoIng and con- 
ventional metal YbCoIng,'"*' in which the superconductivity is 
induced by the multilayers of CeCoIng. The multilayer high- 
Tc cuprates represent a further example of superconductors 
where this type of spin-orbit coupling could play a role . i^-i?) 

We investigated previously the basic properties of multi- 
layer superconductors, such as the electron structure, super- 
conducting order parameters, and magnetic properties with 
focus on the effects of Rashba spin-orbit coupling arising 
from the local non-centrosymmetricity. It has been shown that 
the Rashba spin-orbit coupling significantly affects the super- 
conductivity when the spin-orbit coupling is larger than the in- 
terlayer single particle hopping. Indeed, anomalous behavior 
of the upper critical field has been experimentally observed in 
the superlattice of CeCoIng possibly connected with the para- 
magnetic limiting influenced by spin-orbit coupling.'**' Sev- 
eral other superconductors with locally non-centrosymmetric 
crystal structure have been theoretically investigated. '^"^^' 

In our earlier studies we focused on the superlattice of 
CeCoIng, '''' we neglected the difference of electric potential 



between the inner and outer layers. This approximation may 
be valid for the superlattices of CeCoIng since both inter- 
layer hopping and Rashba spin-orbit coupling can be larger 
than the potential difference. On the other hand, it has been 
shown that the charge imbalance due to the potential dif- 
ference plays an important role for the magnetism and su- 
perconductivity in high-Tc cuprates with more than two lay- 
ers. Thus, it is desirable to clarify the effect of spatially 
inhomogeneous Rashba spin-orbit coupling on multilayer su- 
perconductors with charge imbalance. For this purpose, we 
study the trilayer systems with charge imbalance for simplic- 
ity. 

First, we introduce a model Hamiltonian for two- 
dimensional trilayer superconductors including spatially 
modulated Rashba spin-orbit coupling as 
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where Cksm ('^t„„) is the annihilation (creation) operator for 
an electron with spin s on a layer m. Two outer layers are de- 
noted by m - 1,3, while the inner layer ism -2. We consider 
a simple square lattice and assume the dispersion relation as, 
B(k) - -2f(cos kjc+cos ky)-fi. We choose the unit of energy as 
f = 1 and fix the chemical potential /i = -1, which leads to the 
electron density per site being approximately 0.63 for V = 0. 
The second term describes the Rashba spin-orbit coupling 
arising from the lack of local inversion symmetry. The cou- 
pling constants a,„ should be (01,02,03) - {a,Q,-a) owing 
to the symmetry. We assume a simple Rashba-type form of 
g-vector g(k) = (- sinA;y, sinkx, 0) without going into micro- 
scopic derivation of the antisymmetric spin-orbit coupling. 

The third term introduces intralayer Cooper pairing via an 
off-diagonal mean field. We ignore interlayer pairing as we as- 
sume a weak interlayer coupling. Owing to the spatially mod- 
ulated Rashba spin-orbit coupling arising from the broken lo- 
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cal inversion symmetry, the order parameter Ass'm(k) involves 
both spin singlet and spin triplet components. One of the spin 
singlet and triplet components has the same sign between lay- 
ers, while the other should change the sign. For instance, the 
staggered order parameter of p-wave superconductivity is in- 
duced by the staggered Rashba spin-orbit coupling when the 
d-wave Cooper pairing is predominant as in high-Tc cuprates 
and CeCoIn5. We here neglect this spin triplet component be- 
cause the spin susceptibility discussed below is independent 
of the staggered order parameter as we showed in ref. 13. 
In this paper, we assume the uniform spin singlet order pa- 
rameter with s-wave symmetry {Aii,„(k) = -Aiin,{k) = i/r) 
for simplicity. We have confirmed that the model for the d- 
wave superconductivity gives the qualitatively same results. 
We take {ifrl - 0.01 to be sufficiently small to satisfy the con- 
dition \Ass'm(k)\ « la,,, I <K Bp (sp is the Fermi energy), as 
realized in most (locally) non-centrosymmetric superconduc- 
tors. 

The fourth term in eq. (1) describes the interlayer hopping 
of electrons between nearest-neighbor layers. Since we con- 
sider a quasi-two-dimensional system, we assume an inter- 
layer hopping tj_ = 0.1 is much smaller than the intralayer 
hopping t - 1 . In the following calculations, we will always 
use this value for t_i_, if not stated otherwise. The last term of 
eq. (1) is the static potential V which induces the charge im- 
balance between the layers. We focus on the effect of the po- 
tential V in the following part. Note that the electron density 
per site in the inner layer is «!„ - 0.78,0.68,0.62,0.56, and 
0.49 for V = -0.5,-0.2,0,0.2, and 0.5, respectively, while 
that in the outer layers is Wout ~ 0.63 independent of V. 

We calculate the spin susceptibility for magnetic fields 
along c-axis. The normalized spin susceptibility ;(fs/;tn for var- 
ious potentials V and spin-orbit couplings a are shown in 
Fig. 1, where and;^n are the spin susceptibility at T = 
in the superconducting state and in the normal state, respec- 
tively. The normalized spin susceptibility along a^j-plane is 
only half of that along c-axis as shown in ref. 13. For V = 0, 
our calculation reproduces the results obtained in ref. 13. The 
spin susceptibility vanishes Xs/Xn = for a = as in the 
centrosymmetric spin singlet superconductors. With increas- 
ing a, the spin susceptibility first ascends and shows a peak 
around a ~ tj_, indicating the interplay between spin-orbit 
coupling a and interlayer hopping tj_. Indeed, this is the fin- 
gerprint of crossover from the centrosymmetric superconduc- 
tivity to the non-centrosymmetric superconductivity which 
occurs with increasing the spin-orbit coupling ff.'^' 

When we switch on the potential V, the peak of spin sus- 
ceptibility shifts to a ~ V. This shift does not mean the sup- 
pression of the Rashba spin-orbit coupling. Indeed, for a small 
spin-orbit coupling a < fj_, the spin susceptibility in the su- 
perconducting state is significantly enhanced by the potential 
difference V. This indicates that the effect of Rashba spin- 
orbit coupling is enhanced by the charge imbalance for rea- 
sons discussed below. 

In order to elucidate roles of the charge imbalance in more 
details, we show the contribution of each layer to the spin 
susceptibility in Fig. 2. We assume the positive and negative 
potential V in Figs. 2(a) and 2(b), respectively. First, we dis- 
cuss the small spin-orbit coupling a < tj_- We see that the 
spin susceptibility of outer layers is enhanced by the charge 
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Fig. 1. (Color online) Spin susceptibility in the superconducting state at 
T = as a function of the spin-orbit coupling a. We show the spin suscepti- 
bility along c-axis normalized by that in the normal state. We assume V = 
(solid line), V = 0.2 (dashed-dotted line), V = 0.5 (dotted line), V = 0.7 
(dashed line), and V = 1.0 (dash two-dotted line), respectively. 



the spin susceptibility shows a substantial layer dependence 
for the moderate potential difference V - 0.2, while the spin 
susceptibilities of inner and outer layers are nearly identical 
for y = and a < fj_. This layer dependence appears through 
the decoupling of outer and inner layers due to the charge 
imbalance. As the potential difference \V\ increases, the inter- 
layer hybridization due to tj_ is suppressed. Then, the effect 
of Rashba spin-orbit coupling is enhanced, and therefore, the 
spin susceptibility of the outer layers increases while that of 
inner layer decreases. When the layers are completely decou- 
pled from each other, the spin susceptibility of the outer layers 
is not suppressed by the superconductivity as for uniformly 
non-centrosymmetric superconductors'^^ while that of the in- 
ner layer is completely suppressed as in the centrosymmetric 
superconductor The total spin susceptibility shown in Fig. 1 
is enhanced by the potential difference |y| since the contribu- 
tion from the outer layers is larger than that from the inner 
layer 

Next, we discuss the large spin-orbit coupling a > tj_. We 
see the peak of the spin susceptibility around a ~ \V\m Fig. 1 . 
This peak indicates the crossover in the electronic state. In or- 
der to illustrate this crossover, we show the schematic figures 
of Fermi surfaces for f = (see Fig. 3). The Fermi surface is 
identical for all layers at V = a = as shown by the dashed 
lines in Fig. 3. When the potential of inner layer V > is 
turned on, the Fermi surface of the inner layer shrinks, as 
shown by the solid line in Fig. 3(b). Then, the mismatch of 
Fermi surfaces between m - I and m - 2 and between m - 2 
and m = 3 suppresses the interlayer hybridization. When we 
switch on the Rashba spin-orbit coupling, the Fermi surfaces 
of the outer layers are split. One of the split Fermi surfaces 
having a positive spin helicity is drawn by the solid Unes in 
Figs. 3(a) and 3(c), respectively. The Fermi surface of m = 1 
layer shrinks while that is enlarged for m - 3, because the 
spin-orbit coupling has the opposite sign ai = -03. Then, the 
mismatch of Fermi surfaces between m - \ and m = 2 is 
removed, and therefore, the quasiparticle strongly hybridizes 
between the inner layer m -2 and the outer layer m - I. This 
hybridization enhances the spin susceptibility of inner layer 
without decreasing the spin susceptibility of outer layers, as 
shown in Fig. 2. This is the reason why the total spin sus- 
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Fig. 2. (Color online) Contribution of inner and outer layers to the spin 
susceptibility for (a) positive potential V and (b) negative potential V, respec- 
tively. Thin and thick lines show the spin susceptibility of outer layer and 
inner layer, respectively. We show the results for V = (solid line), |V| = 0.2 
(dashed line), \V\ = 0.5 (dotted line), and |y| = 0.7 (dash-dotted line), respec- 
tively. 
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Fig. 3. (Color online) Schematic figure of Fermi surfaces in the absence of 
interlayer hopping = 0. We show one of the split Fermi surfaces in the 
outer layers. The dashed lines show the Fermi surfaces for a = and V = 0. 
The solid lines show the Fermi surface of m = 1,3 for a > and that of 
m = 2 for V > 0. 



coupling a is furthermore increased, the interlayer hybridiza- 
tion is suppressed again, and then the spin susceptibility of 
the inner layer is decreased. 

In order to substantiate this picture, we now show the single 
particle wave function in the normal state. We focus on the 
single particle state with momentum k = (0, kf) and takes the 
.X-axis for the quantization axis of spin, as in ref. 13. Then, 
the single particle Hamiltonian is block-diagonalized with the 
block 
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the 3x3 matrix, we obtain three bands with the energy 
Ei(k) > E2(k) > E^ik) for the given ^-vector and wave 
function = Zim^;ml ~*)m where /' is the band index. 

Figure 4 shows the representative weights |wj„,p for the band 
7 = 1 on the three layers m. For vanishing spin-orbit coupling 
(or = 0) the single particle weight is localized on one to the 
layers for finite V, i.e. on the inner layer m = 2 in the case of 
band j - 1 . As the spin-orbit coupling is turned on, the lay- 
ers begin to hybridize leading to the strongest hybridization 
for a ~ V (strongest change of weight). Eventually, at larger 
a the weight has shifted to the layer m = 1. Most notably, 
for the strongest hybridization between outer and inner lay- 
ers, the effect of spin-orbit coupling is also transferred to the 
inner layer giving rise to a maximum of the susceptibility. 

Charge imbalance may open a way to observe "local non- 
centrosymmetricity" in multilayer superconductors such as 
high-Tc cuprates even for comparatively weak spin-orbit cou- 
pling by studying the local spin susceptibility. The local spin 
susceptibility is, in principle, accessible in site-selective NMR 
experiment through Knight shift measurements.'^ Experi- 
mental investigations have not focused on the aspect of lo- 
cal non-centrosymmetricity so far. In order to motivate such 
studies, we consider here a simple tight-binding model as an 
approximation to generic band structure of trilayer high-Tc 
cuprates in the overdoped regime in order to avoid complica- 
tions with antiferromagnetic order The inplane hopping in- 
volves nearest- and next-nearest-neighbor hopping, leading 
to s(k) = -2f 1 (cos + cos fcy) + 4f2 cos k^ cos ky - fi and 
the interlayer hopping tj_ is simply nearest neighbor Exper- 
imental studies give for such systems a charge distribution 
among the layers of — 0.9 and Wout — 0.8.1^' These val- 
ues are obtained by the choice of parameters (fi, f2, fx.jf , V) - 
(1,0.3,0.1,-1.08,-0.2) which we now use to calculate the 
spin susceptibility in the superconducting phase for varying 
a as depicted in Fig. 5. We find a clear difference in the spin 
susceptibility between inner and outer layers even for small 
spin-orbit coupling a < 0.04. This is in contrast with the re- 
sult for V - 0. Thus, the moderate charge imbalance in our 
model for multilayer cuprates significantly enhances the ef- 
fect of Rashba spin-orbit coupling on outer layers. Naturally 
the Rashba spin-orbit coupling would also play a role for the 
antiferromagnetic and superconducting states of multilayer 
cuprates. In particular, the difference of magnetic response be- 
tween outer and inner layers would affect the superconducting 
phase structure in the magnetic field. For example, the pair- 
density wave state which has been discussed recently^^' would 
be stabilized by the reduced interlayer hybridization due to 
the cooperation of spin-orbit coupling and charge imbalance. 
A more detailed discussion of the interplay of magnetism and 
superconductivity in the high-Tc cuprates in the light of local 
non-centrosymmetricity will be given elsewhere. 

In summary, we have investigated the eff'ect of charge im- 
balance on the spin susceptibility of trilayer superconductors 
by taking into account the Rashba spin-orbit coupling arising 
from the local non-centrosymmetricity. As an important re- 
sult, we find that charge imbalance enhances the spin suscep- 
tibility of outer layers, because the interlayer hybridization is 
suppressed by the mismatch of Fermi surfaces. The same ef- 
fect leads to a reduction of the spin susceptibility of the inner 
layer. For a small spin-orbit coupling and moderate charge 
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Fig. 4. (Color online) Single particle weight of the band / =1 on the layer (a) m = 1, (b) m = 2, and (c) m = 3, respectively. We assume the potential 

on the inner layer V = (solid lines), V = 0.2 (dash-dotted lines), V = 0.5 (dotted lines), V = 0.7 (dashed lines), and V = 1 (dashed two-dotted lines), 
respectively. 
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Fig. 5. (Color online) Spin susceptibility in the model for the multilayer 
high-Tc cuprates. We show the spin susceptibility of inner (thick dashed line) 
and outer (thick solid line) layer for V = -0.2. The thin lines show the results 
for V = 0. The other parameters are described in the text. 
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layers, thus showing that the effect of spin-orbit coupling is 
enhanced by the charge imbalance. This may also be the case 
for multilayer high- Tc cuprates. When the spin-orbit coupling 
is increased, the strong hybridization occurs between the in- 
ner and outer layers around a ~ \V\, and then, the spin suscep- 
tibility of inner layer increases substantially. These coopera- 
tive and/or competing roles of Rashba spin-orbit coupling and 
charge imbalance can be tested by the experiments for multi- 
layer superconductors such as the high-Tc cuprates as well as 
the superlattice of CeCoIns/YbCoIns. 
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